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ABSTRACT: Staircase and Poisson scaling are two typical patterns we observe for the exclusive 
number of jets at high energy hadron colliders. We examine these scaling properties for photon 
plus jets production at the LHC and find that this channel is well suited to study these features. We 
illustrate and discuss when to expect each of the two patterns, how to induce a transition through 
kinematic cuts, and how photons are related to heavy gauge bosons. Measurements of photon+jets 
production is therefore providing valuable information on exclusive jet scaling, which is going 
to help to eventually understand the theoretical origin of exclusive jet scaling properties in more 
detail. 
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1 Introduction 

After testing and reproducing many interesting aspects of the Standard Model at the LHC, the focus 
of the ATLAS and CMS collaborations is rapidly moving toward searches for Higgs particles [1] 
or physics beyond the Standard Model [2]. The production rates for any of these search channels 
are small, for example compared to W/Z+jets or top pair production, channels which constitute 
their main backgrounds. Traditionally, at hadron colliders we have relied on the appearance of 
leptons, photons or missing transverse energy to point us to interesting new physics processes. In 
this approach QCD effects and jet production are either ignored or considered a nuisance. 

Starting with the suggested searches for a light Higgs boson in weak boson fusion [3], this 
attitude has changed; this search shows how the QCD structure of signal events can be turned into 
a powerful handle to reject large backgrounds. The key analysis tools are (central) jet vetos [4-6], 
which for example suppress QCD-initiated W/Z+jets events or hadronic top pair production. Im- 
plicitly, this approach is adopted in Higgs searches for example in the H — > WW or H — > 77 
channels, when those searches are divided into 0-jet, 1-jet and soon 2-jet strategies [1]. 

In a similar spirit, searches for example for supersymmetry benefit from the measurement of 
the number of jets which includes information on the color structure of the new heavy states [7] ; 
the only caveat is that we need to carefully separate decay jets from QCD jet radiation associated 
with hard processes [7, 8]. What is missing for all such analyses is an experimentally established 
and theoretically sound link between choosing re-jet samples for an analysis and a systematic study 
of the corresponding rtj ets distribution for signal and background processes [5]. A dedicated study 
along this line would map out the behavior of exclusive nj ets distributions after different cuts, under- 
stand its basic features, and quantify the notorious theory uncertainties associated with jet counting. 
As it will turn out, the production of a photon in association with QCD jets is a perfect basis for 
such studies and complements the ones, which have already been performed at the LHC [9]. 

From text book quantum field theory we know that successive p^-ordered photon radiation 
off a hard electron — as well as successive gluon radiation off a hard quark — follows a Poisson 
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pattern for the exclusive number of photons or gluons [10]. This pattern corresponds to a simple 
probabilistic picture of successive independent splitting. The splitting probability is linked to the 
coupling constant, the color factor, the form of the splitting kernel, and a scale logarithm. The 
Poisson scaling pattern, however, has not been observed in inclusive production at hadron colliders 
since UA1/UA2. Instead, for many processes we find staircase scaling, namely a constant ratio of 
exclusive ra-jet rates a n+ i/a n = 12 [11, 12]. This feature has been studied by the LHC experiments 
[13] and also finds application in background-modelling in phenomenological approaches [14]. 

The description of exclusive jet rates is at odds with our description of QCD at hadron col- 
liders. Partem densities obeying the DGLAP equation [15] resum collinear logarithms and absorb 
the corresponding infrared divergences. As a consequence, any computation based on such parton 
densities is jet inclusive, i.e. it allows for an unspecified number of collinear jets radiated off the 
incoming partons. Strictly following the DGLAP approximation, these perfectly collinear jets are 
not observable at the LHC. However, the assumption of perfect collinearity is modified by the ini- 
tial state parton shower, which redistributes the dominantly collinear jet radiation into the physical 
phase space. Evaluating exclusive event samples with exactly n jets is theoretically limited by the 
precision of the parton shower description, including its obvious breakdown for sizeable transverse 
momenta [16]. This is the reason why in the past exclusive nj ets distributions could rarely be ex- 
ploited to compare collider data to QCD predictions. However, as described above, current LHC 
analyses force us to overcome this limitation and study exclusive nj e t s distributions (including cuts 
corresponding to jet veto survival probabilities) starting from perturbative QCD. 

Matching of a hard QCD matrix element with a collinear parton shower [17-20] allows us to 
reliably simulate and study jet scaling patterns from perturbative QCD. It does not require a funda- 
mental re-organization of QCD perturbation theory [21] but simply relies on the proper phase space 
simulation of collinear logarithms. This way it does not only include the radiation of one or two 
hard jets, as correctly described by fixed-order QCD calculations, but any number of radiated jets 
including high multiplicities obviously well described by the parton shower. We use Sherpa [22] 
with its Ckkw [18] matching scheme to describe radiation of up to seven jets with high precision. 
Typically, we check an additional two more jets for unexpected features, but with correspondingly 
reduced statistics. 

As we will show in this paper, a particularly promising channel to measure exclusive jet rates 
and compare them to QCD predictions is hard photon production in association with jets. 

The cross section is large enough to already have enough data to not only validate Monte 
Carlos, but also test potential rej ets scaling hypotheses in various phase space regions. We will 
show how to define and extract different kinematical regimes of the photon to compare various 
hypotheses of QCD radiation with data and to cross check jet emission against other production 
modes. Possible applications towards new physics searches under the premise that Poisson or 
staircase scaling are phenomenologically realized to good approximation have been discussed in 
Refs. [1,4, 5]. 
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2 Staircase scaling 



Some qualitative implications of staircase scaling are already known from the analysis of the inclu- 
sive njets distribution of Standard Model processes. These include pure QCD jets production [7, 23] 
and W/Z production with jets [11, 12, 23]. As we will see, 7+jets production when constrained in 
certain cut scenarios provides an additional channel, hence contributing to a better understanding 
of the possibly observed rij ets scaling patterns. Staircase scaling is defined in terms of constant 
ratios of the experimentally measured jet-inclusive n-jet cross sections 

R(n+i)/n = = R- (jet-inclusive) (2.1) 

It turns out that we can equivalently formulate this condition in terms of inclusive (<j) and exclusive 
(a) numbers of jets. Correspondingly defining 

R{n+i)/n = = R (jet-exclusive) (2.2) 

for the exclusive rates, the resulting ratios are identical [7] , 

K ~ ~~ 1 — ^7 — 1 — — tt ■ U-->) 

°n + L,=n+1 *j °n + <Tn+l L,=0 RJ ^ + R(Jn _L_ 

1 — R 

This way the merits of perturbative QCD and its perturbative predictions directly translate to the jet- 
exclusive final states in a well-defined approach where lower multiplicities are utilized to constrain 
the higher ones. 

While a proper analytical derivation of this feature from first principles is still missing, it is 
observed to a good approximation in both experimental data and theoretical calculations using 
matrix element and parton shower merging [23, 24]. Only using jet radiation via parton showering 
this scaling feature is not matched as well, which is expected given that such inclusive processes do 
not offer a hard scale in relation to which we can define collinear radiation. In addition, for T^+jets 
production it has been shown that fixed order QCD corrections stabilize the observed staircase 
pattern [25]. 

Unlike the other cases mentioned above, at first sight photon production does not posses an 
obvious jet scaling behavior. It only occurs once we include strong separation cuts between the 
photon and each of the jets, effectively removing any logarithmic enhancement linked to QED 
photon radiation. Moreover, as we will show in this paper, different basic cuts can easily induce 
different scaling patterns. 

For our simulation we rely on Sherpa vl.3.0 [22] and its Ckkw matching up to five ma- 
trix element jets. We reconstruct jets using the anti-A^ algorithm from FastJet [26, 27] with 
Raati-k T = 0.4, which gives us a very moderate geometric separation of two jets. When dealing 
with photons in a QCD environment some familiar subtleties have to be considered [28, 29]: a 
photon can arise from non-perturbative fragmentation. Those photons are not useful in our case 
since our focus is obviously not on QED corrections to multi-jet production rates. Therefore, we 
opt for a solid photon isolation. A naive hard cut e.g. on the jet-photon R distance limits the phase 



-3- 



0.5 



0.4 



0.3 



0.2 



0.1 



m, 



m 



i i i 

R„=0.1597 
i yJ >70GeV 

an 

R„=0.1488 
r %j >90GeV 

m=-0.0007 

an 

R„=0.1407 
^j>UOGeV,° 

m =0 .0029 
an 



0.5 



0.4 



0.3 



0.2 



0.1 



Q\ 1 1 1 1 1 

2/1 3/2 4/3 5/4 6/5 7/6 



■ K, 


R„=0.1437 

^-=-0.005 
an 


■ *? 


R.,=0.1266 

" >u -/r 

^■=0.0001 
an 


. R°y 


R„=0.1087 

^-=0.0053 
an 




2/1 3/2 4/3 5/4 6/5 7/6 



Figure 1. Two scenarios establishing staircase scaling for 7+jets production at the 7 TeV LHC. Left: in- 
variant mass criterion of m 7J > 70, 90, 110 GeV for each jet. Right: geometric separation between the 
photon and each of the jets. The extracted values for Rq and dR/dn are defined in Eq.(2.5). The error bars 
correspond to our numerics with 1.6 • 10 7 events. 



space of soft gluon emission and is infrared unsafe. We instead define an isolated photon through 
a hadronic energy deposit of less than 10% of py >7 in a cone of size R < 0.4 centered around the 
photon direction [28]. If this criterion is not met the photon candidate is pushed into the jet finding 
algorithm. 

For reconstructed jets and photons in this section we then require 

p Ta > 50 GeV, p T ,j > 50 GeV, |r/ 7 | < 2.5, \ Vj \ < 4.5 , (2.4) 

where 77 and y denotes the pseudo-rapidity and rapidity, respectively. These cuts are very inclusive 
and democratic, so we can expect to observe the staircase scaling behavior known for pure QCD 
jets. 

We observe that photon plus jets events fulfilling Eq.(2.4) alone do not show any kind of 
simple jet scaling behavior. What we are still missing is the crucial photon-jet separation criterion. 
In Fig. 1 we show two sets of nj e ts distributions for different separation criteria. Inspired by the 
VF/Z+jets analysis we can define a wide photon-jet separation in terms of the invariant mass. In 
Fig. 1 we find that almost prefect staircase scaling appears for minimal values of m 7 j > mz, with 
a very slight degradation for alternative mass scales. The corresponding description in terms of a 
geometric separation leads to very similar results, but only once we require Rjj > 1. This value 
we can understand from the typical m 7J - values combined with px > 50 GeV. In both cases the first 
ratio i? 2 /i is notorious, an effect that has been observed in other channels before [7]. The origin 
lies in a strong PDF suppression of the two-jet configurations. The price to pay for producing an 
additional jet, passing in particular the pxj criterion, is highest for the transition a y j — > cr^jj. 
The measure is the relative increase in partonic center-of-mass energy in order to produce the final 
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Figure 2. Left: effect of a consistent variation of a s (M z ) on staircase scaling. Right: position of the photon 
in comparison to all jets, ordered by px- The two curves are for at least two-jet and at least three-jet events. 

state with one additional jet. This relative increase and the corresponding PDF suppression factor, 
however, becomes rather insignificant for higher jet multiplicities. 

Setting aside the first entry we can test the quality of staircase scaling by fitting the form 

D CTn+l dR 

R(n+i)/n = = #0 + n , (2.5) 

and determine the slope to compare it to the perfect staircase scaling prediction dR/dn = 0. For 
all curves shown in Fig. 1 we find dR/dn in the 0.01 — 0.001 range, essentially compatible with 
zero. The constant values Ro range around 0.14, but with a small spread. In Sec. 4 we will contrast 
these values with the W/Z+jets cases [7, 11, 12]. While the definition of the hard process does 
play a role in determining Rq, in Sec. 3 we will see that the far dominant factor is p™ 1 " fixed in 
Eq.(2.4). 

In the left panel of Fig. 2 we show the effect of varying a s on the observed staircase scaling 
pattern. Between the largest and lowest values of a s there is a 7% difference. The effect of this shift 
on R is correspondingly small. The reason is that a s and for example the gluon parton densities are 
not independently extracted [30]. An increase in the value for the strong coupling is compensated 
by a decrease in the corresponding gluon density, postponing the expected blow up of the theory 
uncertainty to larger n values than we can show in Fig. 2. 

On theoretical grounds, we can link staircase scaling to the presence of the large gluon self 
coupling, i.e. the non-abelian structure of QCD. In the absence of any hard scale from our process, 
relatively hard jets are still mostly generated through initial state radiation (ISR). Our simulation 
confirms that the final state radiation (FSR) cascades initiated by the core process jet and additional 
ISR jets generates the large and democratic jet multiplicities defining staircase scaling. This is 
illustrated in the right panel of Fig. 2. The production of one hard jet with one hard photon is 
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dominated by the partonic subprocess qg qj. Unless we induce a hard scale by cuts, the 
evolution of the incoming quark or gluon is dominated by few and non-abelian splittings after the 
first hard ISR emission. Geometrically, for 2-jet events the quark-photon system typically recoils 
against the harder of the QCD jets. Only from the 3-jet configuration on we can simply split the ISR 
gluon, with a jet separation given by R > 0.4 according to the jet algorithm. Correspondingly, in 
the right panel of Fig. 2 we see that the photon is typically as hard as the hardest jets: if we require 
at least two jets the photon is the hardest or second hardest object in roughly half of the events each 
and the third hardest object only in 0(5%) of all events. For at least three jets in the final state the 
photon is equally likely to be the first, second and third hardest object. In other words, while one 
of the jets might usually recoil against a relatively hard photon, the additional jets responsible for 
the staircase scaling pattern are relatively soft. This is a result of splitting hard ISR gluons. 



3 Poisson scaling 

According to field theory text book knowledge radiating massless gauge bosons off, e.g., a hard 
fermion does not follow a staircase pattern [10]. Instead, successive soft radiation ordered in pt 
yields a Poisson distribution which can potentially be observed in the exclusive number of jets 

e -fi n n 

On = cr : , (3.1) 

n! 

where n is the expected number of emissions. As an example for such a process consider the 
exponentiable purely abelian contributions to multi-jet rates in e + e~ collisions, see e.g. [31, 32]. In 
this section, we construct a cut scenario, for which the approximation of subsequent soft emission 
is sufficiently realized to observe Poisson scaling in 7+jets. 
For the exclusive scaling ratios Eq.(3.1) translates into 

y " a n n+1 

The assumptions entering the derivation of Poisson scaling are twofold: first, there should be one 
splitting function, for example the radiation of a photon or a gluon off a fermion. In the soft 
limit successive gauge boson radiation is automatically ordered by the emission angle. This way, 
we avoid combinatorial factors of the kind n! from differently ordered emission in the numerator. 
The crucial factor 1/n! in Eq.(3.1) appears through the over-counting of the bosonic phase space. 
Poisson scaling is what one expects from a statistical point of view when we assign probabilities 
to statistically independent splittings. An example for this statistical treatment are Sudakov factors 
or collinear splitting probabilities following the DGLAP equation. The reason why solutions of 
the DGLAP equation show a Poisson behavior is that in its derivation we only take into account 
successive splittings of incoming partons on their way from the proton to the hard process. This 
is exactly what corresponds to a resummation of collinear logarithms and the removal of infrared 
divergences through the definition of scale dependent parton densities. The presence of a large 
kinematical logarithm makes parton shower simulations the appropriate tool to reproduce Poisson 
scaling. To force multi-jet scaling for example in association with a photon into such a Poisson 
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Figure 3. Poisson scaling for 7+jets production at the 7 TeV LHC. Left: different transverse momentum 
criteria for the leading jet, all other jets have ptj > 20 GeV. The extracted values for i?o and n are defined 
in Eq.(3.4). The error bars correspond to our numerics with 2.1 ■ 10 6 events. Right: effect of a consistent 
variation of a s (mz). 

regime we can follow this argument of the parton shower and the Sudakov factors. What we need 
is a well defined hard subprocess, e.g. the leading 2-particle 7-ji system which induces many 
successive splittings of the incoming partons. 

Following the distance measures shown in Fig. 1 an obvious choice could be an increased 
value of m 7J1 S> 100 GeV. This is similar to weak boson fusion Higgs production, where the large 
invariant mass of the two tagging jets rrijj > 600 GeV induces a Poisson scaling of the Z+jets 
backgrounds [5] . However, we find that the cleanest Poisson distribution is induced by requiring a 
single hard jet, i.e. requiring 

p r ,7 > 20 GeV, p T ,j > 100, 20, 20, ... GeV, |r/ 7 | < 2.5, \yj\ < 4.5 , (3.3) 

instead of Eq.(2.4). Generating the hard scale through the hardest jet is more efficient than asking 
for a hard photon, because according to our earlier argument a hard photon with pj- a = 100 — 
200 GeV could easily recoil against several jets from splitting ISR. Requiring a hard average px 
for the jets would work as well, though. The important requirement is that through the cuts we 
induce a large scale separation with respect to the pr,min for radiated additional jets. 

In Fig. 3 we see how enforcing a staggered pt cut scenario immediately changes the staircase 
scaling pattern into a Poisson distribution. Already for Pr,n > 100 GeV we see a clear deviation 
from any kind of staircase behavior provided we allow all other jets to be as soft as ptj > 20 GeV. 
For Pr,ji > 150 GeV the leading i? 2 /i rat i° increases to values above unity, which means that in 
the exclusive nj ets distribution the maximum will move away from zero. 

To test the quality of the Poisson description we fit the Rr n+ i> i n distribution which is expected 
to follow Eq.(3.2). If we allow for a deviation from the one -parameter Poisson shape of the kind 

n 

R (n+l)/n = -, + ^0 > (3.4) 
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Rq is reminiscent of the staircase pattern Eq.(2.3) and should come out essentially zero while n is 
the only free parameter in the Poisson shape and gives the expected number of jets. As expected, 
the value of n increases for harder leading jets. Just as for the staircase scaling we do not include 
the first entry i? 2 /i m tne fit- 
While it follows the basic expectation, namely becoming large and exceeding unity, i? 2 /i does 
not fit the Poisson shape well. In the first bin we encounter again a PDF suppression effect for 
producing the second jet. While this effect is much smaller than the mismatch of i? 2 /i f° r staircase 
scaling, it is important to not include this bin into the fit to Eq.(3.4) because it would lead to 
obviously wrong best- fit values for n and Rq. 

If we also exclude the high-multiplicity bins we find small values of \Ro\, closer to zero than 
to i?7/g ~ i?8/7. This is expected. Our argument for Poisson scaling rests on the impact of a large 
scaling logarithm which has to be generated by successive and ordered ISR. Beyond some point 
these successive splittings will stop feeling a large logarithm and we can expect to fall back onto a 
non-negligible staircase scaling. 

For example in the case of Higgs production we know that the large-re limits of R(n+i)/ n i n 
the staircase and Poisson setups show hardly any difference. In Fig. 3 we estimate the staircase tail 
to be in the Rq ~ 0.3 range. This is significantly different from the Rq ~ 0.15 which we find in 
Fig. 1. The reason is simply the reduced px threshold of 20 GeV for the Poisson studies. In the 
right panel of Fig. 3 we again show the consistent variation of a s . Clearly, the dependence is very 
small and does not affect the Poisson scaling feature. 

4 Massive gauge bosons 

Originally, jet scaling studies have been established for iy+jets events. Not surprisingly, Z+jets 
events behave qualitatively and quantitatively the same [7]. The purpose of the first two sections 
of this paper is to show that after strict jet-photon isolation we can see the same scaling patterns in 
photon production. 

Provided that for well separated photons the non-existence of large logarithms between inher- 
ent mass scales leads to staircase scaling we expect 7+jets and Z+jets production to be very similar. 
This very fact is also exploited experimentally to model the invisibly decaying Z background to 
new physics searches [33] from corresponding measurements of 7+jets in control regions [34]. In 
Fig. 4 we first see that based on the same jet cuts as in Eq.(2.4) Z production follows the same stair- 
case pattern. The extracted value Rq = 0. 149 for p-rj > 50 GeV is consistent with the literature [7] 
as well with our findings in Sec. 2. 

For the Poisson regime the situation is slightly different. For the photons a cut on the leading 
jet of 

PT,ji > 100 GeV compared to a reduced general jet threshold of 20 GeV already induces 
a large enough logarithm. The Z mass in the final state could be expected to further enhance this 
scaling logarithm, even though it does not really translate into a collinear logarithm expressed in 
terms of pt j. In the right panel of Fig. 4 we find that the n values we extract from the Z case 
are very close to those for the photons, in particular taking into account the statistical uncertainties 
which affect the result for PT,ji > 200 GeV. The high-multiplicity staircase limit of the Poisson 
distribution for Rq ~ 0.20 — 0.25 again is similar to the photon case. 
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Figure 4. Scaling patterns for Z+jets production at the 7 TeV LHC. Left: staircase scaling for prj > 
50 GcV. Right: Poisson scaling for different transverse momenta of the hardest jet, all other jets have 
Px,j > 20 GeV. The error bars correspond to our numerics with 10 7 events before cuts. 



The agreement of both figures provides a consistency check of our previous statements but 
also opens up the possibility to experimentally cross check Z+jets (and VF+jets) production against 
7+jets in yet another way. This holds not only on the qualitative but also on the quantitative level. 
The only difference between the two channels is that for photons the scaling patterns only appear 
once the photon is very well separated, controlling any additional QED logarithms which for the 
massive Z case do not play any role at these energies. 

5 Outlook 

Counting numbers of exclusive jets has many applications in LHC searches, implemented for ex- 
ample as distinct Higgs analyses for different rij ets values [1] or central jet vetos [3]. To apply 
such cuts while maintaining stability of the theoretical predictions we need to properly understand 
exclusive nj ets distributions both experimentally and theoretically. A major obstacle to overcome 
is that exclusive jet measurements have no straightforward interpretation in fixed order perturba- 
tion theory and exclusive quantities are typically plagued with large theoretical uncertainties. On 
the other hand, multi-jet merging, for example using the Ckkw scheme [18], allows us to study 
njets distributions including a free choice of kinematic cuts and has proven successful in various 
experimental analyses so far. 

Unlike VF/Z+jets and pure QCD jets production, the associated production of jets with a hard 
photon naively does not show simple scaling patterns. We show that once we require a widely 
separated photon we recover staircase scaling cr n+ i /cr n = Ro for the total cross section [7, 11, 12]. 

Once we induce a large logarithm through kinematic cuts we see how the scaling pattern turns 
into a Poisson distribution for the exclusive number of jets. This is known for weak-boson-fusion 
cuts in Higgs production [5]. For our photon channel a transverse momentum cut on the leading jet 
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(and not on the photon) works best. Only in the high-multiplicity regime an underlying staircase 
pattern remains. 

Given our observations and the large available photon sample at LHC, this channel is espe- 
cially well suited to study jet scaling, including a proper experimental and theoretical error analysis. 
Combined with more channels where we expect staircase scaling, such a comprehensive study will 
not only provide crucial information on MC validation but will also help to eventually reveal the 
origin of staircase scaling. 

A translation of these photon measurements into W/Z+jets production is straightforward and 
will significantly impact new physics searches: typically, 7+jets production is used to infer (Z — > 
z^)+jets production in both control and signal regions [34]. Therefore, a better understanding 
of exclusive jet quantities and the translation of 7+jets into (Z — > z^z/)+jets production can help 
to systematically reduce the uncertainty of the background extrapolation in e.g. the jets+missing 
energy channel. 
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